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ABSTRACT 



Current-driven reconnection is investigated as a possible mechanism for an- 
gular momentum transport in astrophysical disks. A theoretical and computa- 
tional study of angular momentum transport from current-driven magnetohydro- 
dynamic instabilities is performed. It is found that both a single resistive tearing 
instability and an ideal instability can transport momentum in the presence of 
azimuthal Keplerian flow. The structure of the Maxwell stress is examined for a 
single mode through analytic quasilinear theory and computation. Full nonlinear 
multiple mode computation shows that a global Maxwell stress causes significant 
momentum transport. 



1. Introduction 

As matter accretes from astrophysical disks onto a central body (such as protostars, 
neutron stars, and black holes) angular momentum is rapidly transported outward. This 
redistribution of angular momentum is far too rapid to be explained by collisional viscosity. 
It is believed that turbulence initiated and sustained by magnetohydrodynamic (MHD) 
plasma instabilities can be responsible for the enhancement of effective viscosity in accretion 
disks (Shakura & Sunyaev 1973; Pringle 1981). Disks appear to be hydrodynamically 
linearly stable according to the Rayleigh criterion. However, it is not clear whether through 
a nonlinear process a hydrodjTiamically stable flow can become turbulent, a topic that has 
been debated vigorously. Recently, it has been shown experimentally that hydrodynamic 
turbulence in astrophysically relevant flows at Reynolds numbers up to several million 
cannot provide the required effective viscosity (Ji, Burin, Schartman & Goodman 2006). 

The ineffectiveness of hydrodynamic turbulence in transporting momentum has 



- 3 - 



motivated the search for other instabihties that could yield enhanced viscosity. A candidate 
for such instability is the magnetorotational instability (MRl). The MRI, originally studied 
by Velikhov 1959 and Chandrasekhar 1960, was analyzed for astrophysical disks by Balbus 
and Hawley in the early 1990s (Balbus & Hawley 1991). They showed that a differentially 
rotating disk in the presence of weak magnetic field is MHD unstable. Thus, under 
appropriate conditions, the MRI can lead to turbulent momentum transport in disks (see, 
e. g., Balbus & Hawley 1998, and the references therein). The momentum transport arises 
from Maxwell or Reynolds fluid stresses. In this paper, we investigate whether those stresses 
and momentum transport can be provided by a current-driven MHD instability, rather 
than through the flow-driven MRI. In particular, we consider in depth a tearing instability, 
which is a resistive MHD instability, that causes magnetic reconnection and can persist in a 
strong magnetic field (in regions in which the MRI could be stabilized by the strong field). 

Magnetic fields have long been considered to be prevalent in some disk regions, as 
well as to play a key role in jets emanating from disk regions. Observational evidence is 
in hand for the existence of strong magnetic fields in disks, although the observations are 
limited. In recent years, direct measurements have shown strong evidence of magnetic field 
in protostellar disks. Existence of a significant azimuthal magnetic field of about 1 kG has 
been confirmed by observations around FU Orionis in the innermost regions of its accretion 
disk (Donati et al. 2005). The complex topology of magnetic field around stars during their 
formation has also be confirmed through observation and extrapolation of the reconstructed 
surface magnetic field. Recent observation and analysis of magnetic field on the surface of 
classical T Tauri stars (cTTS) suggest strong octopolar field 1.2 kG) and smaller dipolar 
field (^ 0.35 kG) (Donati et al. 2007). 

There is some indirect evidence of strong magnetic field in protosolar nebula based 
on the magnetization of primitive meteorites. It was suggested that an MHD nebular 
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dynamo can generate a magnetic field with intensity of up to lOG at tlie distance of a 
few astronomical units from the Sun (Levy 1978), and was shown that this magnetic field 
can produce flares in a nebular corona (Levy & Araki 1989), which may account for the 
magnetized meteorites. 

The source of large-scale magnetic fields of cTTS could be the combination of fossil field 
(from interstellar media) and field generated through dynamo action (Donati et al. 2007). 
However, the source of large-scale magnetic field in disks is still unknown. Large-scale 
magnetic fields in disks may be supplied externally and may also be generated through a 
kinematic alpha-omega dynamo (Reyes-Ruiz & Stepinski 1999; Torkelssson & Brandenburg 
1994; Rudiger, Elstner & Stepinski 1995). MRI can also produce a turbulent small-scale 
MHD dynamo (Brandenburg et al. 1995; Hawley et al. 2001), however, the generation 
of large-scale magnetic field by MRI is not known to occur. Based on the observational 
evidence of magnetic field and the modeling of self-generated large-scale magnetic field 
(MHD dynamo) in disks, we examine current-driven instabilities as a candidate for angular 
momentum transport in an MRI stable magnetized disk with large currents. 

In this paper, we explore the possibility of angular momentum transport in disks by 
current-driven instabilities in situations where MRI is stable (due to strong magnetic field) 
and can not provide the required effective viscosity. We recognize three classes of situations 
that current-driven instabilities may be important in transporting angular momentum in 
astrophysical disks; 1) in the inner disk region around a young magnetized protostar 2) in 
the innermost and active regions of weakly-ionized protostellar and protoplanetary disks 3) 
in the upper and lower surfaces layers of a protoplanetary disk and magnetized coronas. 
Below we further discuss these situations. 

The strong and complex large-scale magnetic fields detected in protostellar disks can 
provide the free energy (current) for current-driven instability and magnetic reconnection 
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in these systems. Disk-star interaction and accreting matter to the stars along the complex 
magnetic field lines have been studied through MHD simulations by Long, Romanova 
& Lovelace (2008) (and the reference therein). However, with open and complex field 
lines attached to the disk observed at radii only up to several stellar radii and with disk 
differential rotation, magnetic reconnection may be important in the inner region of a disk 
close to a young protostar (such as cTTs) and cause radial angular momentum transport. 

Current-driven instabilities may also be important in protoplanetary and protostellar 
disks (e.g. protosolar nebulas, disks around CTTs). However, in weakly ionized 
protoplanetary disks, the magnetic field might not be well coupled to the gas throughout 
the whole disk. Therefore, magnetic reconnection can only play a role in some parts of these 
disks which are thermally ionized (the innermost region r<0.1AU) or nonthermally ionized 
by cosmic rays (in the outer active region). Other possible candidates for momentum 
transport in disks are dynamo modes and fiow-driven MRL The notion of angular 
momentum transport through magnetic stresses from the oscillating dynamo modes has 
been examined by Stepinski & Levy (1990) in a protostellar nebula in the absence of 
fiow-driven MRI turbulence. It was estimated that a magnetic field strength of 10^ — lO'^G 
is needed for transport of angular momentum from protostar into a surrounding disk. 
Torkelssson & Brandenburg (1994) also found stronger magnetic torques from quadropolar 
dynamo solutions. The effectiveness of MRI in transporting angular momentum transport 
in weakly-ionized disks is limited and strongly depends on the ionization factor. It has been 
shown that in the so-called dead zone between the active layers in protoplanetary disks, 
MRI can not operate (Gammie 1996; Fleming & Stone 2003; Zhu, Hartmann & Gammie 
2010). However, including nonideal effects such as Hall effect and ambipolar diffusion have 
shown to change the stability and the saturation of MRI in weakly-ionized protoplanetary 
disks (Wardle 1999; Balbus & Terquem 2001; Wardle 2007). 
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In strongly magnetized coronas where MRI is stable, the role of current-driven 
instabilities can become significant. Magnetic activity at the surface of the disks, possibly 
originated from the internal dynamos, may produce coronas, which can transport angular 
momentum vertically and radially. Using local simulations. Miller & Stone (2000) showed 
the generation of strongly magnetized corona (force-free) through an MRI-driven turbulence 
originated in the core of a disk with initially weak magnetic field. They also investigated 
angular momentum transport from MRI in a vertically extended local domain. However, 
the possibility of angular momentum transport in the strongly magnetized, MRI stable 
force-free corona was not studied. Theoretical models have also been developed to study the 
possibility of angular momentum transport by the coronal magnetic field (Goodman 2003; 
Uzdensky & Goodman 2008). The current-driven instabilities can also become important 
in corona, upper and lower layers of protoplanetary disks, with significant currents, where 
MRI neither operates in the core region (dead zone) nor is unstable in the coronas. 

In this paper, we specify the equilibrium magnetic field (and corresponding current 
density) in disk geometry, which is chosen to be sufficiently strong that the fiow-driven 
MRI is stable. We then study current-driven instabilities. Unlike fiow-driven MRI and 
pressure-driven interchange instabilities, which can be treated locally, current-driven 
instabilities are driven due to equilibrium gradients with scale length comparable to the 
global scale lengths. It is shown that gradients in global equilibrium parallel current provide 
the driving source for both ideal and resistive current- driven (tearing) instabilities. Ideal 
current-driven instabilities are global modes (with global velocity perturbations) with fast 
Alfvenic growth rates which are independent of resistivity. However, for tearing modes, 
resistivity is important in the vicinity of the radius at which the wave number parallel to 
the magnetic field vanishes (i.e., /c|| = k ■ B/B = 0). Magnetic reconnection (or tearing) 
occurs at such radial locations. Although the reconnection occurs locally in radius, the 
tearing modes are global in their radial extent. Therefore, tearing modes are usually treated 
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using asymptotic matching where solutions outside the reconnecting layer (outer layer or 
ideal 77 = solutions) are matched with the inner layer (reconnecting layer rj ^ 0) solutions. 

We investigate the physics of momentum transport through three sets of calculations of 
increasing completeness. First, we perform linear stability analysis to establish conditions 
for instability and to investigate radial structure of the modes. Second, through quasilinear 
analytic theory, we calculate the Maxwell stress in the outer ideal region and show its 
dependence on the azimuthal flow shear and instability growth rate. We find that in the 
presence of azimuthal Keplerian flow, the radial and azimuthal magnetic fluctuations are in 
phase and causes a nonzero Maxwell stress for a single (a single spatial Fourier component) 
current-driven mode and thus momentum transport. We also examine the structure of the 
stress for a single mode in the quasilinear regime. We show that the stress and the resulting 
momentum transport is localized to the vicinity of the reconnection layer. However, an 
ideal current- driven instability has a global Maxwell stress, and therefore produces more 
global transport. 

Third, we compute the simultaneous nonlinear evolution of multiple modes, each 
corresponding to different reconnection layers (separated radially). With the inclusion of 
nonlinear mode coupling, the Maxwell stress becomes global, extending over the full plasma 
cross-section. We perform nonlinear computations which are 1) weakly nonlinear driven 
2) strongly nonlinear driven. In the weakly nonlinear driven case, many current-driven 
modes are linearly unstable while the effect of nonlinear coupling to stable modes is weak. 
In the strongly nonlinear case however, many stable modes are nonlinearly driven by 
linearly unstable modes. It is shown that in the weakly nonlinear driven case, the nonlinear 
structure of Maxwell stress is global mainly due to the growth of an ideal current- driven 
instability. However, in the strongly nonlinear case, the global structure of the Maxwell 
stress is mainly due to the nonlinear coupling of many tearing modes. We also estimate 
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the effectiveness of the current-driven instabihties for angular momentum transport using 
the standard Shakura-Sunyaev a model for both the weakly nonlinear driven and strongly 
nonlinear driven cases. We find that Maxwell stress is much larger for the strongly nonlinear 
case compared to the weakly nonlinear case and ass are about 10"'^ and 10"^ for the weakly 
nonlinear driven and strongly nonlinear driven case, respectively. 

The paper is organized as follows. The characteristics of momentum transport 
by current-driven instabilities and the model used are described in sections |2] and |3l 
respectively. Single mode calculations are presented in section SJ In Section 14. 1^ equilibrium 
and linear stability analysis are described. Analytical quasilinear calculation of Maxwell 
stress and the structure of single mode in quasilinear regime are presented in section 14.21 
We present full nonlinear multiple mode computations in section |5l Multiple computations 
for a weakly nonlinear driven and a strongly nonlinear driven case are presented in sections 
15.11 and |5.2[ respectively. We summarize in section |6l 



2. Physical attributes of momentum transport by 

current-driven instabilities 

The momentum transport caused by current-driven instability has the same origin 
in fluctuation-induced stresses as does the well-studied flow-driven MRI. In that sense, 
the mechanism for transport is the same, although the underlying MHD instability is 
different in its energy source and spatial structure. The change in flow is given by the MHD 
momentum equation 

P^ = -pV.VV + JxB (1) 

For a rotating, cylindrical plasma we decompose the flow and magnetic fields into spatially 
mean and fluctuating quantities, where mean values (denoted by <>) are averaged over 
axial and azimuthal directions and fluctuating values subtract out the mean. If we then 
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average the above equation over the axial and azimuthal directions we find 

d < \ > 

p = -p< V.VV > + < J X B > (2) 

at 

which describes the time evolution of the mean fiow which depends upon radius and time. 
We see that the mean fiow can evolve from fiuctuating fiow (arising from a Reynolds stress, 
the first term on the RHS) and fiuctuating magnetic field (the Maxwell stress or Lorentz 
force of the second term on the RHS). These terms, quadratic in the fiuctuations, can 
be nonzero upon spatial averaging. MHD current-driven instabilities produce spatially 
fiuctuating quantities that, through the Reynolds and Maxwell stress alter the radial profile 
of the axial and azimuthal mean fiow. These forces can decrease the fiow at one radius, 
while increasing it another, essentially transporting momentum radially. In this sense, the 
mechanism is identical to momentum transport from the well-studied magneto- rotational 
instability. Crucial differences lie in the energy source, spatial structure, and dynamics of 
the underlying instabilities. 

The MRI is a global ideal instability (occurs at zero resistivity) with the stresses acting 
globally in radius. The current-driven tearing instability also has an amplitude that is 
global in radius, but the forces on the RHS of the Eq. [2] act locally in radius for a single 
instability. For a single tearing instability magnetic reconnection occurs at one radial 
location that for which the wavelength parallel to the equilibrium magnetic field becomes 
infinite (Biskamp 1993). That is, the instability amplitude is constant along the helical 
magnetic field at this radius, although it varies in the azimuthal and axial directions. Only 
the vicinity of this radius is electrical resistivity important. The resistive layer of narrow 
radial extent, within which reconnection occurs, functions as a boundary layer between two 
ideal regions. Within this region amplitudes of velocity and current density can become 
large, leading to locally strong forces that alter fiow. 

An additional difference between the MRI and tearing instability is that momentum 
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transport in the MRI can be understood as the instabihty acting to diminish its energy 
source (flow gradient). For the tearing instabihty the analogue is the reduction of the 
current density, which is a strong effect of the instabihty. Reduction of the flow gradient 
can be viewed as a parasitic effect of the tearing mode. Interestingly, the MRI also alters 
the magnetic field (Ebrahimi, Prager & Schnack 2009), so both instabilities alter the flow 
and field. We should note that an ideal current-driven instability despite having the same 
energy source as tearing instability, exhibits different characteristics. A major difference 
is that for an ideal current- driven instability (with no reconnecting layer in the plasma 
volume), the fluctuation- induced forces are global as will be shown below. 

The tearing mode, being a resistive instability, does not lend itself to simple analytic 
calculation of instability-induced transport. However, in section I4l2l an analytic calculation 
of the transport is provided for the ideal region. We find that in the presence of azimuthal 
flow, the radial and azimuthal magnetic fluctuations are in phase and results in a non-zero 
Maxwell stress. The resulting stress (and the direction of angular momentum transport) 
depends on the global equilibrium as well as the global radial mode structure. We therefore 
employ quasilinear computations to show the structure of stresses and outward momentum 
transport. It will be shown that the structure of the stresses for ideal current-driven 
and tearing modes are very different, therefore they can affect the momentum transport 
differently. With multiple tearing modes present, nonlinear three- wave interactions alter 
the structure of the modes and the forces broaden radially. 

3. The model 

Throughout this work, we employ the MHD equations in doubly periodic (r, 0, z) 
cylindrical geometry for both the analytic and computational studies. All variables are 
decomposed as f{r,(j),z,t) = /m,fc('^, ^)e*^""^+''^^ =< f{r,t) > +f{r,(f)z,t), where 
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< / > is the mean [m = k = 0) component, and / is the fluctuating component (i.e., 
all other Fourier components with m 7^ and k ^ 0). Note that the mean component 
is a function of radius. We consider an azimuthal equilibrium flow Vq = V^{r)(j) in a 
current-carrying disk configuration plasma. (Here, V^{r) is the equilibrium value of the 
mean azimuthal flow < >). In order to excite the current-driven instabilities, both 
vertical magnetic field and the azimuthal magnetic field, B = Bz{r)z + B^{r)(j), are imposed 
(Fig. [1]). The single fluid MHD equations are, 
(9A 

_ = -E = 5VxB-r/J (3) 
= -SpV.VY + S3xB + P„y^Y - S^VP - 5pV$ (4) 

CM 2 



dP 

'dt 
dp 

'dt 

B = V X A (7) 



-SV ■ (PV) - S{T - 1)PV ■ V (5) 
-SV ■ (pV) (6) 



J = V X B (8) 

where the variables, p, P,V, B, J,, F, and $ are the density, pressure, velocity, magnetic 
field, current, ratio of the specific heats, and gravitational potential respectively. Time 
and radius are normalized to the resistive diffusion time = 47rr2^/c^?7o and the outer 
radius r2, making the normalized outer radius unity. Velocity V and magnetic field B are 
normalized to the Alfven velocity Va, and the magnetic field on axis Bq, respectively. The 
parameter 5* = — is the Lundquist number (where ta = ?"2/Va), and P^ = measures 
the ratio of characteristic viscosity uq to resistivity 770 (the magnetic Prandtl number), 
where Tms = 47rr2^/c^i^o is the viscous diffusion time. The factor Pq = SttPq/Bq is the beta 
normalized to the axis value. The resistivity profile rj is uniform. The boundary conditions 
in the radial direction are as is appropriate to dissipative MHD with a perfectly conducting 
boundary: the tangential electric field, the normal component of the magnetic field, and the 
normal component of the velocity vanish, and the tangential component of the velocity is 
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the rotational velocity of the wall. The azimuthal (0) and axial (2) directions are periodic. 

We pose an initial value problem that consists of the equilibrium plus a perturbation 
of the form f{r,(j),z,t) = fm,k{r,t)exp{i'm(f) + ikz). Equations (1-6) are then integrated 
forward in time using the DEBS code. The DEBS code uses a finite difference method 
with a staggered grid for radial discretization and pseudospectral method for azimuthal 
and vertical coordinates. In this decomposition, each mode satisfies a separate equation of 
the form dfm,k/dt = Lm,kfm,k + Zl(m',fc') ^m,k,m' ,k' ^ where L^.fc is a linear operator that 
depends on fQfi{r,t), and N^j^^' ^' is a nonlinear term that represents the coupling of the 
mode (m, k) to all other modes (m', k'). ( This latter term is evaluated pseudospectrally.) 
The time advance is a combination of the leapfrog and semi-implicit methods (Schnack et 
al. 1987). 

We present three types of initial value computations, linear, nonlinear single mode 
and fully nonlinear computations. In linear computations, the initial conditions consist 
of an equilibrium < /(r) > plus a single mode /m,fc(?", O)exp(2m0 + ikz) perturbations. 
Only the mode (m, k) is then evolved; in particular, the equilibrium (the m = 0, k = 
mode, or < / >) is not evolved, and remains fixed in time. In the nonlinear single mode 
computations, however, the m = 0, k = component (the background) is allowed to evolve 
self-consistently. The evolution of the background profile /o,o can affect the evolution of 
the mode (m, k) and cause the mode to saturates. In a fully nonlinear computation, all 
modes are initialized with small random amplitude and are evolved in time, including the 
full nonlinear term (iV^^fc^^'^^')- 

The plasma rotates azimuthally with a mean Keplerian flow < Vtp{r) >oc r"^/^. 
The initial (i.e., at t = 0) radial equilibrium force balance (Eq. 2) is satisfied by 
^Vp + pV$ = pV^/r, where V$ = GM/r'^, and a magnetic force free condition J x B = 0. 
The initial pressure and density profiles are assumed to be radially uniform. Pressure and 
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density are evolved; however, they remain fairly uniform during the computations. 

We consider a cylindrical disk-shaped plasma with aspect ratio L/(r2 — ri) (where L 
is the vertical height, and ri and r2 are the inner and outer radii, respectively) (Fig. [1]). 
The inner and outer radial boundaries are perfectly conducting, concentric cylinders that 
can rotate independently at specified rates. Periodic boundary conditions are used in the 
vertical and azimuthal directions. 

The aspect ratio (L/(r2 — ri)) used in the nonlinear computations is 1.3. The nonlinear 
computations are performed in a thick-disk approximation where vertical and radial 
distances are of the same order. The range of parameters used in the computations is 
/3 = 1 - 10, 5 = 10^ P,„ = 0.1 - 20. 

4. Single mode calculations 

Here we present linear and quasilinear single mode computations for current-driven 
instabilities. We first choose equilibria which are unstable for current-driven instabilities and 
then perform linear computations (section 14.11) . To study whether a single current-driven 
instability transports momentum, we examine the Maxwell stress in the ideal MHD region 
(outer region) through quasilinear analytical calculations (section l4.2p . We also examine 
the structure of quasilinear stresses for a single mode (with specific m and n) with the 
equilibrium described in section 14.11 

Unlike flow-driven MRI and pressure-driven interchange instabilities, which can 
be treated locally, current- driven instabilities have global characteristics. Gradients in 
global equilibrium parallel current provide the driving source for the ideal and resistive 
current- driven instabilities. Ideal current-driven instabilities, so-called kink modes, are 
helical long wavelength structures with global (broad radial extent) velocity perturbations. 
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These ideal modes do not scale with resistivity and have fast Alfvenic growth rates. In an 
ideal plasma, the restoring force is infinite and fluid is frozen to the field. However, in the 
presence of resistivity, the field lines can break up and at locations where the parallel wave 
number is zero, (so called resistive or reconnecting layer) field lines can reconnect. Both 
structure and growth rate of a resistive current-driven mode (tearing mode) are affected 
by resistivity. The solutions around the reconnecting layer are solved by including the 
resistivity (so called inner layer solutions) and are asymptotically matched with solutions 
outside the reconnecting layer (so called outer layer solutions). We will show that the 
structure of the stresses for ideal current-driven and tearing modes are very different, 
therefore they can affect the momentum transport differently. 

4.1. Linear computations 

We first choose equilibrium profiles that are unstable for current- driven tearing 
instabilities. We start with a force free plasma J x B = 0. Current flows parallel to the 
magnetic field line J|| = A(r)B, where J\\ is the parallel current. Figure |2]^a) shows a 
typical equilibrium A(r) = J||/B profile used in the computations at t=0. The equilibrium 
vertical and azimuthal magnetic field components can be obtained from the force balance 
equation V x B = A(r)B, and are shown in Fig. [2]J^a). The magnitude of the magnetic fields 
is large enough to make this equilibrium MRI stable. 

As described in the previous section, perturbations are assumed in the form of 
exp(ikr) ~ exp(im0 + i27rnz/L) in cylindrical geometry , where m and n are the azimuthal 
and axial mode number respectively [kz = 27m /L). At locations where parallel wave 
number is zero k ■ B = mB^/r + 27rnBz/L = (so called resonant or reconnecting surfaces), 
resistivity becomes important and reconnection can occur. From this resonant condition, a 
field line winding number q = —m/n = 27TrBz/LB^ can be defined. The profile of axial 
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winding number q is shown in Fig. [2l Here, conventionally positive q corresponds to a 
negative axial mode (n) and vice versa. It will be shown that depending on the equilibrium 
profile (A), current-driven modes with —5 < q < 5 can be linearly unstable. 

To verify that the equilibrium is tearing mode unstable, we first perform linear single 
mode computations (single m and n) in the absence of azimuthal flow. The radial mesh 
points used is 250. The modal structure and the growth rates of the tearing modes depends 
on the equilibrium profile, i.e parallel current profile A(r). Therefore, we perform four set 
of linear computations with four different equilibria. The four equilibria are distinguished 
by the maximum value of the equilibrium A = J\\/B profile. We find that both resistive 
current- driven modes (tearing modes) and ideal current- driven modes with azimuthal 
modes numbers m=0,l,2,3 are linearly unstable for sets of equilibria with A^ax =10.6 
(shown in Fig. [2]), Amax =14.2, and Amax =17.2. For equilibria with Amax =9, only m=0,l,2 
tearing modes are linearly unstable. Below, we further discuss the mode structure and the 
properties of tearing modes and ideal current- driven modes. 

In the presence of resistivity, current- driven instabilities (tearing) can become linearly 
unstable. The radial structure of radial magnetic and velocity eigenf unctions for m=l 
tearing mode with axial modes number (n=-l) are shown in Fig. [3]^a) (for A^ax = 10.6) and 
Fig. [3](b) (for Amax = 17.2). As it is seen, the reconnecting component of magnetic field 



Table 1: Growth rates {'JTa) of tearing modes (m,n), for two A, 



Amax 


(1,0) 


(1,1) 


(1,-1) 


(0,1) 


(2,1) 


(2,-1) 


(3,1) 


(3,-1) 


9 


stable 


0.030 


0.031 


0.0362 


0.0155 


0.0165 


stable 


stable 


10.63 


0.058 


0.059 


0.045 


0.057 


0.048 


0.0445 


0.029 


stable 


14.2 


0.23 


0.148 


0.101 


0.132 


0.141 


0.0557 


0.116 


0.0105 


17.2 


0.3 


0.31 


0.21 


0.275 


0.3 


0.14 


0.255 


0.082 
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(Br) is nonzero around the reconnecting surface [r=0.6 for (l,-l)in Fig. [3]|^a)]. Moreover, 
the radial fluctuating velocity is concentrated and change sign around the reconnecting 
surface. One of the characteristics of tearing instability is the jump in the logarithmic 
derivative of Br across the resistive layer (A' = (-B^lrs ~ ^r\rs)/^r\rs), where rs denotes 
reconnecting surface.. For tearing mode to be unstable, this jump should be positive 
(Furth et al. 1963). The growth rate of the tearing mode scales with resistivity and A' as 
7tcaring oc r/'^/^A'^/^ (or oc S~^^^A'^^^). As is seen in Fig. [3l A' is positive for m=l tearing 
modes, and is larger in Fig. |3t^b) with larger free energy (current gradient Amax), which 
results in larger growth rate. We also note that for larger Amax? (Irl) mode becomes a 
double tearing mode (reconnects at two locations), and radial velocity perturbation changes 
sign twice in radius (absolute values of the eigenfunctions are shown in Fig. [3](b)). 

The growth rates of m=0-3 tearing modes for four set of A profiles with Amax=9, 10.6, 
14.2, and 17.2 are given in Table 1. For smaller Amax? the tearing mode growth rates are 
smaller, and the growth rates increase with Amax approaching the ideally unstable limit 
(approaching Alfvenic growth rates). Modes with azimuthal mode numbers 0-3 and axial 
mode numbers 1,-1 have tearing mode structure with a corresponding reconnecting surface. 
However, mode with axial mode number zero (1,0), is ideally unstable for equilibria with 
Amax=10.6, 14.2, and 17.2. The radial structure of this mode is shown in Fig. H] which has 
a kink-like characteristic. Unlike the tearing modes, ideal current-driven modes does not 
necessarily have to have a reconnecting surface (with k ■ B) within the plasma and the 
fluctuating velocity components can be very global without changing sign (kink-like). The 
growth rates of the ideal current-driven modes also don't scale with resistivity. It has been 
shown that the current-driven modes are unstable for the equilibria chosen (Fig. [2]) and 
therefore, we can study the effect of both tearing and ideal current-driven modes on the 
momentum transport through nonlinear computations (section [5]) . 
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4.2. Analytical quasilinear calculations 

In order to investigate momentum transport from current-driven instabilities, we obtain 
MHD stresses using quasilinear calculations of a single mode with an initial azimuthal 
Keplerian flow. The question is whether a single current-driven mode can transport 
momentum and affect the azimuthal flow profile. To obtain more insight into this question, 
we first analytically examine the ideal MHD equations with azimuthal flow in cylindrical 
geometry. We then construct the Maxwell stress < r'^B^B^ >= 2Re{r'^B^Br) from the 
linearized solutions in the outer (ideal) region. Through analytical quasilinear calculations, 
we aim to identify whether the Maxwell stress is nonzero for a current-driven instability 
and therefore momentum can be transported. We do not intend to analytically solve the 
complete sets of solutions (inner and outer solutions), and we only present the simplified 
outer solutions. To construct the radial structure of the quasilinear stresses, we use the 
eigenf unctions from the linearized computations. 

The linearized incompressible ideal MHD equations in the presence of mean azimuthal 
flow are, 

+ (Vo ■ V)V + (V ■ V) Vo) = -V(Bo ■ B) + (Bo • V)B + (B • V)Bo (9) 

_ = V X (V X Bo) + V X (Vo X B) (10) 

We assume perturbation of the form Q(r, 0,z,t)=Q(r) exp(i w t + i (m + /c z)), 
and the equilibrium magnetic field and flow are B = B^{r)e^ + Bz{r)ez and Vo = V(j,{r)e^ 
respectively. The linearized equations can be combined and be presented in the form of one 
ordinary differential equation for \4, 

(rK)" + ai(rV;)' + a2(rK) = (11) 



- 18 - 



where, 



1 dA mG 2mnujVs 
al = ; ttttt + 



A dr ru{l - M^) r^F^{l - M^) 

1 ,d , mGA , d , ImB^F 
a2 = — — 



and; 



A dr ru{l — M^) dr ruj{m? + /c^r^) 
rco r ra;(l — M^)(m^ + /c^r^) ru 
^^2^2(^2 + p^2)^W ^£l;2(m2 + Pr2)(l-M2)V 

^ F(m2 + A;2r2) (1 - M"^)' r oo{m'^ + k^r'^){l - A/Py 



^ = ( 2^.2 2V G = {1- M')V; - (1 + M')V^/r 



w = a; + mV'^/r, M = y/p^ 

r 



F = + = -mB,/r + kB^ 



where u = Uj. + i'j, and Ur and 7 are the real frequency and the the growth rate, 
respectively. The equation for parallel magnetic field perturbations = -^11 terms of 
radial velocity perturbation is obtained. 



fill - ^s[^.(rV,) + ^^^^ _ ^^^^^ )(rV.)] (13) 

To obtain the Maxwell stress term, we write the parallel magnetic field perturbation 
Eq. [13] in terms of radial magnetic field perturbation Br = FVr/oJ and the parallel magnetic 
field perturbation without mean flow B^^\ 

B„ ^ (1 - MYnr - a,H^^^^^)irS,.) + .M^,_J^^2L_„.S,) ,14) 



- 19 - 



where sf°^= i[Fr(r5^)' - r{k±B)\{rBr)]/iB {rn? + klr'^)) is the parallel magnetic field 
in the absence of mean fiow obtained from EqJT3| and also known from Newcomb equation 
(Newcomb 1960), and \ = J\\/ B . In the absence of mean fiow, it can be seen that B^^'^ is 
purely imaginary and B^ and B^^^ are out of phase, and therefore parallel and azimuthal 
Maxwell stresses vanish for the ideal region. However, with the azimuthal mean fiow the 
Maxwell stress < r'^B^B.,. >= 2Re{r'^B^Bj.) from a current-driven mode is nonzero and is 
obtained, 

< r%B,. >= 2, (2*pr^(i;4jMi|s(»)|B, 

+2, ( 4kpr \ '^-^TftfL ^fl.g; ) (15) 



-27 (2kpr ^^^^2 ^ k^r^)FB^ 



where, 



kBu iB, d ~ , 

The azimuthal Maxwell stress Eq. [T5] has been calculated using Eqs. [T31 and Uni assuming 
Br is purely real. With mean azimuthal fiow, Eq. [15] shows that the joint effect of mode 
growth and mean flow produces a nonzero stress term in the outer region. We note that 
for an ideal current-driven mode without a reconnecting surface (Ft^ everywhere), the 
stress term can be global. The sign and the structure of the resulting stress depend 
on the global equilibrium as well as the global radial mode structure and do not follow 
from Eq. [15] alone. However, for a special case of an ideal current-driven mode m=l, 
n=0 (k=n/R=0), discussed in section [^?T| we can examine the sign of the Maxwell 
stress using local WKB approximation. For this mode the Maxwell stress is simplifled as 
< r'^B^Br >= 2Re{[^-^{rBr)]*Br}, where Br is complex, and we have used Eq. [TB] for 
the toroidal magnetic perturbation (or equivalently from V ■ B = 0). Using a local 
approximation, a negative Maxwell stress < r'^B^Br >= —2r^^Bl is obtained, where kr 
is the local radial wavenumber. The negative Maxwell stress results in a bi-directional 
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Lorentz force and causes outward momentum transport. Although WKB is not a vahd 
approximation for the current-driven instabihty and does not provide the structure of 
the transport, it demonstrates an outward transport, which is consistent with the global 
solutions given below. 

To obtain the radial structure of stresses, we employ quasilinear computations. Using 
the linear eigenfunctions for Br (shown in Fig. H]), the radial structure of the Maxwell 
stress term for an ideal current- driven mode, m=l,n=0, is shown in Fig. |5l Similarly, the 
quasilinear Reynolds stress term, < r^V^K >, can also be constructed from the linearized 
solutions for the m=l, n=0 mode (shown in Fig. Oby dashed-dotted line). We note that 
the total stress term, < r'^B^Br > — < r^V^V^ >, is negative everywhere, which causes a 
transport of angular momentum outward. Total azimuthal fluid force consisting of Lorentz 
and inertia terms < J x B >0 -p < V ■ VV >^= < r'^B^B,. > -pj^ < r'^V^Vr >], is 
bi-directional and causes transport of momentum outward, as shown in Fig. [51 

For a resistive current-driven mode (a tearing mode) since the growth rate scales 
as S~^^^ or r/^/^, where S is the Lundquist number, the Maxwell stress term is small in 
the outer region (77 = 0) and the main contribution arises from the inner layer solution 
(resistive region). Equation [15] only presents the outer solution and is singular around the 
reconnection layer (F=0). To obtain the radial structure of stresses for a tearing mode 
with a reconnecting surface k ■ B = F = 0, the inner layer equations need to be solved 
(Ebrahimi, Mirnov, & Prager 2008). We therefore perform nonlinear single tearing mode 
computations with two set of equilibria (A = 10 and A = 17). The radial structure of the 
quasilinear Maxwell and Reynolds stress terms during the linear phase for two cases are 
shown in Fig. [6] The mode structure for the case with Xmax = 10 is a single mode tearing 
structure. As discussed, the main contribution of MHD stresses for the tearing mode arises 
in the inner region, which is also confirmed by the computations [Fig. [6t^a)]. As is seen 
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in Figure ini^a), the MHD stresses are localized around the reconnecting surface (r ^ 0.6), 
and are small in the outer ideal region. Because the computation is in the visco-resistive 
regime, the Maxwell stress is much larger than the Reynolds stress. The localization of the 
stresses leads to a localized flattening of the azimuthal flow within the resistive layer. It 
is also seen in Fig. that the azimuthal Lorentz force {< J x B >0) is bi-directional 
and transports momentum outward. As the free energy for the instability increases, for 
the larger current [Xmax = 17) case, double tearing mode become unstable which has two 
reconnecting surfaces at two radii (r ^ 0.5, r ^ 0.8). The radial structure of the Maxwell 
stress shown in Fig. MJo) is broader than the single tearing mode case. 

It should be mentioned that a more generalized set of ideal compressible equations 
including equilibrium flows in a cylindrical current-carrying plasma (Bondeson, lacono 
& Bhattacharjee 1987) and the gravitational force (Keppens, Casse & Goedbloed 2002; 
Blokland, et al. 2005) have been obtained. They used Frieman & Rotenberg (1960) 
formalism to obtain first order differential equations for displacement and 11 (total kinetic 
and magnetic pressure perturbations). Using these generalized equations (Eqs. 14-20 in 
Blokland et al. 2005), after some algebra we have also obtained the parallel magnetic 
field perturbation (B ■ B/i?) which reduces to Eq. [13] in the incompressible limit. We 
therefore, find that in the incompressible regime by ignoring the pressure perturbation (and 
with uniform equilibrium pressure), but including the gravitational force perturbation, the 
resulting Maxwell stress term remains the same as Eq. [151 

We also note that in a current-free plasma with a uniform axial magnetic field and 
a differential rotation, Eqs. [11] and [12] reduce to a differential equation for a global 
axisymmetric m=0 MRI (Velikhov 1959 and Chandrasekhar 1960), 



+ -IT - - + ^ + T-2 — ^ - r^^^ = 17 

ar^ r or (o;^ — 7^) (07^ — 7"''^ 



where, Q = V^{r)/r, B = BqC^, and u\ = k B^/p. It can be shown that the stresses from 
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an axisjTumetric MRI can be global and depend on flow shear (Ebrahimi, Prager & Schanck 
2009). 

5. Multiple mode computations 

Here, we investigate full nonlinear dynamics of momentum transport from current- 
driven instabilities in a disk geometry using multiple tearing mode computations. In the 
nonlinear computations, with multiple modes (in both azimuthal and vertical directions), 
the additional effect of nonlinear mode coupling is revealed. Moreover, the transfer of 
energy from fluctuations to the mean field can occur during the nonlinear saturation. 
We perform full nonlinear computations when the current-driven modes with different 
azimuthal and axial mode numbers are included in the computations. Current-driven modes 
can nonlinearly interact and cause nonlinear growth. Three sets of nonlinear computations 
are performed in a disk geometry. The first two multiple mode computations are weakly 
nonlinear driven and discussed in section 15.11 A multiple mode computation, which is 
strongly nonlinear driven, is presented in section 15.21 In the weakly nonlinear driven case, 
many current-driven modes are linearly unstable while the effect of nonlinear coupling 
to stable modes is weak. In the strongly nonlinear case however, many stable modes are 
nonlinearly driven by linearly unstable modes. 

5.1. Weakly nonlinear driven 

Here, two sets of computations are performed in a disk geometry with similar force-free 
equilibrium shown in Fig. [2] with (1) Xmax = 9 and (2) Xmax = 14.2. Both computations 
start with a Keplerian flow profile with an on axis amplitude of V^/Va = 0.8 (with Pm = 1, 
S = 10^, f3o = 10, and radial, azimuthal and axial resolutions nr=220, < m < 21 and 
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—43 <n < 43, respectively). 

The computations start with a current-carrying equihbrium, and the free energy from 
the parallel current causes the current-driven instabilities to grow. The radial magnetic 
energy for different current-driven tearing modes are shown in Fig. [71 Tearing modes which 
are linearly unstable (with the growth rates given in Table 1 for equilibria with Xmax = 9 
and Xmax = 14.2), are also driven linearly in the nonlinear computations. Moreover, other 
modes with higher azimuthal mode numbers are driven nonlinearly. As seen in Fig. El^a), 
modes with tearing parity m=0,l,2 (both n=l and n=-l) linearly start to grow and saturate 
around t=0.015 r^j. For tearing mode m=l,n=l, there is a second nonlinear growth before 
saturation. Modes with higher azimuthal mode numbers also shown, and is seen that m=3 
mode linearly grow with small growth rate. This mode is linearly stable in the absence 
of flow and here with the azimuthal flow becomes linearly unstable. For the case with 
equihbrium Xmax = 9, the current gradient (free energy) is not large enough to make the 
ideal current-driven modes (including non-resonant) unstable. Thus for this equilibrium, 
only resonant resistive tearing modes are unstable. 

When the free energy, i.e current gradient, increases both resistive and ideal current- 
driven modes becomes unstable. Figure [7] (b) shows the radial magnetic energies for 
several current- driven modes with the equilibrium Xmax = 14.2. As seen, tearing modes 
m=0,l,2,and 3 linearly grow as expected from linear stability analysis (Table 1). These 
tearing modes also have a nonlinear growth around t= 0.002r/j. In addition to linearly 
unstable modes, m=4 tearing mode which is linearly stable, nonlinearly starts to grow 
around t=0.002r/j. For this equilibrium, as shown in Fig[71^b) ideal current-driven mode 
m=l ,n=0 mode is linearly unstable and saturate about two order of magnitude higher 
amplitude. 

As current- driven modes grow and nonlinearly saturate through modifying the source 
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of instability, i.e the current gradient, they also modify the mean azimuthal flow profile. 
Figure [H] shows the modification of the azimuthal flow profile during the nonlinear evolution 
for both cases Xmax = 9 and Xmax = 14.2. For the case with equilibrium Xmax = 9, the 
modification of flow profile is mainly due to the resistive tearing modes. The azimuthal 
flow profile at two times ti = 0.012, before saturation, and t2 = 0.023, after the saturation 
of tearing modes, are shown in Fig. IHl^a) with Xmax = 9. At t=ti, the flow is modified 
around r=0.6 due to m=l tearing mode. As discussed before, the azimuthal flow is modified 
through the Maxwell and Reynolds stresses. The structure of the fluid stresses are shown in 
Figure ini The Maxwell stress at t=ti is localized around r=0.6 and cause the modification 
of the flow and transport of momentum outward. After the nonlinear saturation, the 
Maxwell stress from all the tearing modes becomes broader (Fig. [U^b)), thus cause a broader 
flow modification as shown at t=t2 Fig. IHl^a). The structure of total nonlinear Maxwell 
stress (from multiple tearing computation) is broader than the structure for a single tearing 
mode (Fig. |9|^a)). The structure is more concentrated in the plasma core (around r=0.5 to 
r=0.8), where the current gradient is large. 

The momentum transport from current-driven modes becomes stronger when ideal 
modes are also present. The nonlinear evolution of azimuthal flow profile for the case 
with larger current gradient X^ax = 14.2 is shown in Fig. MJo). The structure of Maxwell 
and Reynolds stress terms are also shown in Fig. [TOl As seen, the Maxwell stress is 
concentrated in the core region at t=ti before the nonlinear saturation (Fig fTOT a)) and the 
flow modification also occurs around the region where the stresses are peaked from r=0.4 to 
r=0.7. Around the nonlinear saturation at t=t2 and after the saturation at t=t3, it is shown 
that fluid stresses become more global (Fig. [Tor b).(c)) and cause a global modification 
of the azimuthal flow profile and momentum transport outward (Fig. [8]^b)). The global 
modification is due to both resistive tearing and ideal current- driven modes. 
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It is also worth estimating the effectiveness of the current- driven instabihties for 
angular momentum transport using the standard Shakura-Sunyaev a model. In this model 
the effective viscosity is parametrized by the quantity ass which relates to the Reynolds 
and Maxwell stress terms by < pVrV^ > — < BrB^ >= ass < P >, where P is the 
average pressure. The normalized alpha used for our simulations is ass/l^of^- For the two 
weakly nonlinear driven cases presented here, Amax =9 and Aniax=14.2, the time-averaged 
Shakura-Sunyaev a, < ass >, during nonlinear saturation are < ass >~ 5 ■ 10~^ 
and < ass >~ 4 ■ 10""^, respectively. We have also calculated, the time-averaged 
Shakura-Sunyaev a for a weakly nonlinear case with Amax=14.2 and /3o = 1, which is 
< ass >~ 2.5- 10~3_ 



5.2. Strongly nonlinear driven 

We further investigate the effect of nonlinear mode coupling in a different equilibrium 
setting. The computations are performed in a force-free equilibrium for which the 
components of mean magnetic fields are shown in Fig. [11] (solid lines). Full nonlinear 
computation starts with a Keplerian flow proflle with an on axis amplitude of V^/Va = 0.8 
(with Pm = 5, S = 10^, P = 10, and radial, azimuthal and axial resolutions ^^=250, 
< m < 11 and —43 < n < 43, respectively). The free energy, parallel current, is mainly 
concentrated in the inner half plasma region and causes a current-driven instability with 
mode number m=l, n=-l to become linearly unstable. The radial magnetic energy for 
different current-driven tearing modes is shown in Fig. [121 Almost all the modes are 
stable up to around t/r^j = 0.004 and plasma is in a quasi-single mode state (with m=l, 
n=-l). Two other tearing modes, m=l, n=-2 and m=2, n=-l also grow linearly with 
small growth rates. As seen in Fig. [T2](a). around t/r^j = 0.004, due to nonlinear mode 
coupling axisymmetric m=0 mode becomes nonlinearly unstable and saturates at a large 
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amplitude comparable to the amplitude of the initial m=l linearly-driven mode. Other 
non-axisymmetric modes are also driven nonlinearly and a turbulent state is formed. The 
magnetic spectrum for m=0-2 during the nonlinear state is shown in Fig. [T^ b). as can be 
seen, the spectrum shows a broad range of magnetic fluctuations. The transition from a 
quasi-single mode to a multiple-mode state occurs due to strong mode-mode coupling and 
transfer of energy from m=l modes to m=0 modes. In the weakly nonlinear driven case 
(section l5.ip . the equilibrium was linearly unstable for most of the current-driven modes, 
and the nonlinearly driven modes (m=3 and m=4 modes for the cases with Xmax = 9. and 
^max = 14.2, respectively) would not grow to large amphtudes. Here most of the tearing 
modes are driven nonlinearly and saturate at large amplitudes. 

Tearing fluctuations during the nonlinear state (t/r^ > 0.005) can affect the mean 
profiles through the fluctuation-induced convolutions terms. Mean magnetic fields during 
the nonlinear saturation, which are affected by the magnetic fluctuations, are shown in 
Fig. [11] (the dashed lines). As it is seen, the vertical magnetic field changes sign and 
toroidal flux is redistributed by the fluctuation-induced term, < V x B >, the so called 
dynamo term. The tearing fluctuations also affect the flow profile and cause momentum 
transport. The radial structures of Maxwell and Reynolds stresses during the two states, 
quasi-single mode and multiple modes, are shown in Fig [T31 During the quasi-single mode 
state [t/rji < 0.004), the Maxwell stress transports momentum outward, but it is localized 
in the inner region around the reconnection layer. However, during the multiple-mode 
state (t/r^? > 0.005) the structure of total stresses is very global and all the tearing modes 
contribute to the momentum transport. It can also be seen that the Maxwell stress is much 
stronger (about an order of magnitude) for the strongly nonlinear driven case compared to 
the weakly nonlinear driven case (Fig. [TOj) . 

We have also performed strongly nonlinear driven cases with the same force-free 
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equilibrium shown in Fig. [TT]but with Pm=20, /3 = 10, and Pm=20, /3 = 1. The structure 
of stresses after nonhnear saturation (after the growth and saturation of nonhnearly driven 
mode, m=0) are shown in Fig. (TH As seen. Maxwell stress is broad due to nonlinear mode 
coupling. We have calculated the time-averaged Shakura-Sunyaev a during the nonlinear 
saturation, which is < ass >~ 4.3 ■ 10"'^ and < ass >~ 3.3 ■ 10"^ for the two cases with 
/5 = 10 and /3 = 1, respectively. 



6. Summary 

We have shown that tearing instabilities can transport angular momentum effectively 
in thick disks. A single tearing mode transports momentum locally, in the vicinity of the 
radial reconnection location. However, with multiple modes simultaneously unstable and 
nonlinearly coupled to each other, the resultant Maxwell stress and momentum transport 
can span the entire plasma cross-section. It is also shown that an ideal current-driven 
instability without a reconnecting surface has a global Maxwell stress and therefore can 
produce more global momentum transport. However, the global Maxwell stress from 
many tearing modes causes much stronger momentum transport compared to an ideal 
current- driven instability with a global stress. These conclusions are drawn from a hierarchy 
of calculations, including linear stability, quasilinear analytic calculations, and computation 
of the full nonlinear evolution of multiple modes. 

These instabilities, although present in rotating plasmas, are driven magnetically by 
the plasma current. We have selected for study disk equilibria for which the magnetic 
field is sufficiently strong that the plasma is stable to the flow-driven magnetorotational 
instability. Determination of the relevance of tearing instabilities to astrophysical disks 
requires extension of these studies to a wider range of conditions, such as to thin disks and 
different magnetic field strength. As the disk becomes thinner, the relevant tearing modes 
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(with smaller wavelengths) may become more stable. The effectiveness of tearing modes in 
transporting momentum in disks also depends on the ratio of magnetic energy to the flow 
energy. Although in the linear regime, ideal current-driven instabilities are independent 
of resistivity, nonlinear momentum transport from both tearing and ideal current-driven 
instabilities in the low magnetic diffusivity regime could be an interesting topic for a future 
work. 

The authors wish to thank useful discussions with Dalton Schnack, Ellen Zweibel and 
Amitava Bhattacharjee. 
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Fig. 2. — Radial profiles of (a) equilibrium azimuthal and vertical magnetic fields (b) A 
and axial winding number q= 27rrBz/ LB^. 
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Fig. 3. — Radial structure of linear magnetic and velocity eigenfunctions for m=l unstable 
tearing modes, S = W, Pm = ^ (a) Xmax =10.6 (b) Xmax =17.2. 
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Fig. 4. — Radial structure of linear magnetic and velocity eigenfunctions for the ideal m=l, 
n=0 current-driven mode. 
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Fig. 5. — Radial structure of quasilinear total stress (solid line), Maxwell stress < r'^B^Bj. > 
(dashed line), Reynolds stress — < r'^V(j,Vr > (dashed-dotted), and total azimuthal fluid 
force < J X B >0 —p < V ■ VV >^ (dashed-dotted-dotted) constructed from the linear 
eigenf unctions for ideal current-driven mode m=l n=0 mode. 
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Fig. 6. — Radial structure of Maxwell, Reynolds stresses and the Lorentz force during the 
linear phase before the mode saturation for m=l, n=-l (a) single tearing mode with Xmax = 



10 and (b) a double tearing mode with Xmax = 17, S = 10^, Pm = 1- 
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Magnetic energy Wm,n = 1/2/ Br(m,nydr^ vs time for different tearing modes 



(m,n) for (a) Xmax = 9 and (b) Xmax = 14.2 
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Fig. 8. — .The evolution of flow profile during nonlinear computation for (a) Xmax = 9 at 
two times ti=0.012, t2=0.023 and (b) Xmax = 14.2 at three times ti=0.0028, t2=0.0038, and 
t3=0.0056. 
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Fig. 9. — Radial profiles of total (surface-averaged including all the modes) Maxwell and 
Reynolds stress terms at two times (a) ti=0.012, (b) t2=0.023 {Xmax=Q) 
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Fig. 10. — Radial profiles of total Maxwell (solid lines) and Reynolds stress (dashed lines) 
terms at three times (a) ti=0.0028, (b) t2=0.0038, and (c) t3=0.0056, A„ax=14.2. 
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Fig. 11. — Equilibrium azimuthal and vertical magnetic field profiles for strongly nonlinear 
driven case at t=0 (solid lines) and at t=0.0062 during nonlinear saturation (dashed lines). 
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Fig. 12. — (a) Magnetic energy Wm,n = 1/2 / Br(m,n)^dr^ vs time for different tearing modes 
(m,n), (b) magnetic energy spectrum; for the strongly nonlinear driven case. 
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Fig. 13. — Radial profiles of total Maxwell (solid lines) and Reynolds stress (dashed lines) 
terms (a) at ti = 0.0047 during the single mode state (b) at t2 = 0.0062 during the nonlinear 
multiple mode state. 




Fig. 14. — Radial profiles of total Maxwell and Reynolds stress terms after nonlinear satu- 
ration for (a) /3o = 10 (b) /3o =1; Pm =20, S=10^ 



